A Classical Cepheid in a LMC eclipsing binary: evidence of shortcomings in 
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ABSTRACT 

The recent discovery and analysis of a classical Cepheid in the well detached, double- 
lined, eclipsing binary OGLE-LMC-CEP0227, has provided the first determination of 
the dynamical mass of a classical Cepheid variable to an unprecedented 1% accuracy. 
We show here that modern stellar evolution models widely employed to study Galactic 
and extragalactic stellar systems, are able to match simultaneously mass and radius 
(and effective temperature) of the two components with a single value for the age of 
the system, without any specific fine-tuning, assuming the typical metallicity of LMC 
Cepheids. Our conclusion is that there is no obvious discrepancy between dynamical 
and evolutionary masses for the Cepheid star in this system, contrary to previous claims 
of an overestimate of the Cepheid mass by stellar evolution theory. 

Subject headings: binaries: eclipsing — stars: evolution — stars: fundamental parame- 
ters — stars: variables: Cepheids 



1. Introduction 

Classical Cepheid variables have been the subject of several theoretical and observational stud- 
ies, due to their fundamental role in establishing the extragalactic distance scale (see, e.g. Friedman 
et al. 2001; Sandage et al. 2006, and reference therein). The empirical determination of their pul- 
sational properties provides also sets of observational benchmarks for testing both pulsational and 
stellar evolution models. One well known problem related to the theory of classical Cepheids is 
the discrepancy between mass estimates obtained from pulsational and stellar evolution models. 
Specifically, by comparing results of theoretical pulsation calculations with the period, the esti- 
mated mean absolute magnitude and effective temperature of a Cepheid variable, a 'pulsational 
mass' can be derived. At the same time, a comparison of stellar evolution models with just the 
estimated mean absolute magnitude and effective temperature, a so-called 'evolutionary mass' can 
be assigned to the same object. Early investigations by Stobie (1969), Cogan (1970) and Rodgers 
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(1970) discovered that these two independent mass estimates did not agree, the pulsational masses 
being smaller by 20-40%. This disagreement has been named the 'Cepheid Mass Discrepancy'. 

The first fundamental step towards solving this problem was discussed by Andreasen (1988), 
who suggested that an increase of a factor of ~ 2.5 of the radiative opacities in the temperature 
regime 1.5 < (T/IO^ K) < 8 would be able to remove almost completely the discrepancy. This 
increase was nicely confirmed years later by accurate opacity computations performed by the OPAL 
group (Rogers &: Iglesias 1992). Following the implementation of these new opacities in stellar 
model computations, there has been a large number of studies devoted to re-investigate the mass 
discrepancy problem (see, e.g. Keller &: Wood 2002, 2006; Caputo et al. 2005; Evans et al. 2007, 
and references therein). These analyses concluded that the discrepancy was significantly reduced 
to about 10-15% in mass. 

Attempts to further reduce the disagreement from the point of view of stellar evolution mod- 
elling have invoked the inclusion of mass loss, and an increase of the convcctive core mass during 
the core H-burning stage (see, e.g., Caputo et al. 2005; Bono et al. 2006; Keller 2008, and references 
therein) compared to the values determined from the Schwarzschild criterion. Keller (2008), and 
more recently Hilding et al. (2010), have shown that mass loss does not appear a viable solution 
because the required mass loss efficiency would not only be extremely large, but also its depen- 
dence on the stellar mass would be at odds with independent empirical constraints. Hilding et al. 
(2010) have also presented compelling evidence to support the idea that a moderate increase of 
the convective core size during the main sequence would be able to eliminate completely the mass 
discrepancy. The proposed increase of the convective core mass may be achieved as a consequence 
of a change in the adopted physical inputs, and/or the inclusion of additional physical processes like 
rotational mixing and mild convective core overshooting. Major modifications to the stellar model 
input physics do not seem realistic (see, e.g. Cassisi 2010), whereas a host of observational data 
- independent of Cepheid stars - on, e.g., eclipsing binary systems with at least one component 
on the main sequence or near the turn off, the turn off morphology in color-magnitude-diagrams 
of open clusters, star counts along the main sequence and central He-burning sequences in star 
clusters (see, e.g. Andersen 1991; Pietrinferni et al. 2004, and references therein), demonstrate the 
need to extend the main sequence convective cores beyond the formal Schwarzschild boundary. The 
added extension of the convective core size, A„.o! is often parametrized in terms of Acco = ^Hp, 
where Hp is the pressure scale height at the Schwarzschild boundary and A is a free parameter. 
The current, more widely used libraries of stellar models all include an extension of the convective 
core size beyond the Schwarzschild boundary (see, e.g., Pietrinferni et al. 2004; Bertelli et al. 2008; 
Dotter etal. 2008, and references therein). 

Pietrzyhski et al. (2010) have very recently provided the first accurate determination of the 
dynamical mass of a classical Cepheid variable in a well detached, double-lined, eclipsing binary in 
the Large Magellanic Cloud (OGLE-LMC-CEP0227). Thanks to the availability of a high-quality 

photometric dataset, the spectroscopic follow-up, and the near-perfect characteristics of the system 
to derive accurate masses for both components, Pietrzyhski et al. (2010) were able to estimate the 
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mass of the pulsator to an unprecedented 1% precision. The discovery and analysis of this binary 
system is a pivotal step forward in testing the accuracy of both pulsation and stellar evolution 
theory. Pietrzyhski et al. (2010) emphasize that their dynamical estimate of the Cepheid mass 
is in good agreement with the estimated pulsational mass, whilst they suggest the existence of a 
discrepancy with the mass predicted by stellar evolution models. In view of the high precision of the 
parameters of this binary system, and the huge effort devoted in this last decade towards improving 
stellar evolution models, we consider that the estimate of the Cepheid evolutionary mass for this 
system deserves a more specific investigation. The accurate determination of masses and radii for 
both components allows us to test for the first time in a more direct way the accuracy of theoretical 
evolutionary masses of Cepheid stars. Following a standard and powerful methodology of eclipsing- 
binary studies (i.e., Andersen 1991), we will test whether in the age-radius diagram, evolutionary 
models calculated for the empirical values of the mass of the two components and a typical LMC 
composition, can match the observed radii for a common value of the age. With this procedure we 
want to address the following question: is it possible to reproduce the mass/radius constraints of 
both components of OGLE-LMC-CEP0227 by employing a set of widely used off-the-shelf modern 
stellar evolution models, without any fine-tuning for this specific object? 

The next section describe briefiy the models employed in our analysis and presents a comparison 
with the results for OGLE-LMC-CEP0227. A final discussion closes the paper. 

2. The comparison between stellar evolution theory and observations 

We have employed the same state-of-the-art input physics and stellar evolution code adopted 
for building the widely used BaSTI stellar model archive^ (see Pietrinferni et al. 2004, for details), 
to compute a series of models with masses between ~4.1 and ~4.2 Mq, appropriate for the OGLE- 
LMC-CEP0227 components. The BaSTI archive has been extensively tested against observations of 
local stellar populations and eclipsing binary systems (see, e.g. Tomasella et al. 2008, and references 
therein), and is widely employed in studies of Galactic and extragalactic resolved star clusters (see, 
e.g. De Angeh et al. 2005; Mackey & Broby Nielsen 2007; Marfn-Franch et al. 2009, for just a 
few examples), in the determination of Star formation history and chemical enrichment history 
of resolved Local Group galaxies (see, e.g., Carrera et al. 2008), and also in studies of integrated 
properties of extragalactic stellar populations (see, e.g.. Carter et al. 2009). 

We adopted in our calculations the typical initial chemical composition of LMC stars, Z=0.008 
and Y=0.256 - that corresponds to [Fe/H]=— 0.35 for the solar (Z/X) ratio by Grevesse & Noels 
(1993) - e.g. the same composition used by Pietrzyhski et al. (2010) to determine the pulsational 
mass of the Cepheid variable. The convective cores along the main sequence have been extended by 
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an amount Agco = 0.2Hp, as in the BaSTI models for this mass range. As discussed in (Pietrinferni 
et al. 2004) the cahbration of Acco in the BaSTI models is based on the morphology of color- 
magnitude-diagrams of open clusters, and is completely independent of the evolution in the Ccphcid 
instability strip. Finally, all models presented in this work have been computed by adopting a 
Reimers (1975) mass-loss law, with the free parameter rj set to 0.4, and a solar mixing length 
calibration, as in the BaSTI database (Pietrinferni et al. 2004). 

Our analysis detailed below shows that the BaSTI stellar evolutionary models widely employed 
to study Galactic and extragalactic stellar populations are able to match the mass and radius of 
both components of the binary system studied by Pietrzyhski et al. (2010) with a single value for 
the age of the system, without any specific fine-tuning. This implies that there is no significant 
discrepancy between dynamical and evolutionary masses for the Cepheid star in this system. This 
is a result opposite to the conclusions by Pietrzyhski et al. (2010), who suggest an overestimate of 
Cepheid masses by stellar evolution theory. 

The relevant parameters of the eclipsing binary system derived by Pietrzyhski et al. (2010) 
are reported in table 1. The values of the stellar radius and T^ff for the Cepheid (the primary 
component A, as labelled in Pietrzyhski et al. 2010) are mean values suitable to be compared with 
results from stellar models in hydrostatic equilibrium. 

Figure 1 displays a comparison theory-observations in the age-radius diagram. This diagram 
ensures the most direct test between theoretical and empirical quantities, without the need to apply 
color - Tgff relations and bolometric corrections to the output of stellar evolution calculations. 

A satisfactory fit with stellar evolution models of the appropriate mass must reproduce - within 

the error bars - the radii of the two components for a common age. The three selected values for 
the intial mass of the models (Mj) and the corresponding mass during the central He-burning 
phase (Mf) are all within the la error bar of the empirical values. Mass loss implemented with 
the Reimers formula is efficient mainly along the red giant phase before central He-ignition, and 
reduces the stellar mass by a constant amount, equal to ^O.OSM©. Analogous computations with 
77=0 (i.e. with inefficient mass loss) provide essentially the same results in the R — age diagram, 
the only difference being that in this case Mj = Mf by definition. 

Even a cursory analysis of Fig. 1 shows that stellar evolution models provide masses in complete 
agreement with the results from the eclipsing binary analysis. Considering the model with Mj =4. 14 
Mq (M/=4.11 Mq) for the Cepheid component, and the one with M,=4.2l Mq (M/=4.18 Mq) 
for the companion, the observed radii are reproduced within the la error bars for a common age 
of ~154 Myr. These results demonstrate that there is no discrepancy between evolutionary and 
dynamical mass of the Cepheid star in this binary system, when stellar evolution models include 
updated input physics and the extension of the convective cores satisfy independent constraints 
from open cluster color-magnitude-diagrams. 

Figure 2 displays an additional comparison in the R — T^ff diagram. The stellar models 
displayed are the same as in Fig. 1 and the agreement with observations is again excellent for 
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the same combinations of Mi values that satisfies the R — age diagram. We have plotted also the 
Fundamental Blue edge (FBE) and Fundamental Red edge (FRE) of the Cepheid instability strip as 
provided by Fiorcntino ct al. (2007) pulsational models, for a metallicity of Z=0.008, and a mixing 
length equal to 1.5Hp, their reference choice for the treatment of the superadiabatic layers. The 
two components appear to have been captured - as also remarked by Pietrzyhski et al. (2010) - in a 
very short-lived evolutionary stage, i.e. when they move along the brighter branch of the blue loop 
during the core He-burning stage. By examining Fig. 2 we can add that the A component is located 
at the boundary between the fundamental and the first overtone region within the instability strip, 
whereas the B component appears to have left the instability strip since a short time. The fact 
that the B component is not a variable star, and its location with respect to the FRE, poses an 
important constraint on the efficiency of superadiabatic convection in the pulsational models: any 
substantial decrease of the mixing length adopted in the pulsational models would shift the FRE 
towards lower Tg.ff (e.g. Fiorentino et al. 2007), moving the B component within the instability 
strip, at odds with the observations. 

Before closing this section we wish to comment briefly about the effect of our assumptions about 
metallicity and mass-loss, for there are no spectroscopic estimates of the chemical composition of 
OGLE-LMC-CEP0227 components, and no strong constraints on the mass-loss rate of Cepheids 
and their progenitors. Our analysis, as well as the value determined by Pietrzyhski et al. (2010) 
for the pulsational mass of the Cepheid component, is based on the assumption of a 'typical' iron 
abundance [Fe/H]=— 0.35, that is consistent with the mean value determined spectroscopically by 
Romaniello et al. (2008) for a sample of 22 LMC Cepheids. These authors also find a spread 
of [Fe/H] values for their sample, with [Fe/H] ranging between —0.62 and —0.10 dex. Given 
the lack of [Fe/H] estimates for OGLE-LMC-CEP0227 components, we have tested to what extent 
assumptions on their metallicity are critical for the outcome of this analysis. We found that varying 
[Fe/H] by ±0.10 dex (i.e. varying Z between 0.006 and 0.01, and Y using the He-enrichment ratio 
AY/AZ=1.4, as in the BaSTI archive) around the reference value still ensures consistency between 
theory and observations in the age-radius diagram. Values of [Fe/H] outside this range prevent 
models - with mass consistent with the empirical values - from crossing the Cepheid instability 
strip during the central He-burning phase. As a consequence, one would need in this case to fine- 
tune the evolutionary calculations to match this system, for example by including overshooting from 
the bottom of the convective envelopes, that favours the development of loops in the Hertzsprung- 
Russell diagram during the central He-burning phase. 

As for the mass-loss assumptions, we have already remarked that neglecting mass loss or em- 
ploying the Reimers prescription with 77=0.4 (a standard choice for the value of this free parameter) 
do provide consistency between evolutionary and empirical masses. Observations have not deter- 
mined yet conclusive and stringent bounds on mass-loss rates for Cepheids and their progenitors. 
Very recently, Marengo ct al. (2010) have found some indications that the Cepheid prototype 6 
Cephei may be currently losing mass, at a rate in the range of ~ 5 10^^ to 6 1O~^M0 yr^. As 
an extreme test, we have calculated models employing the upper limit of this mass-loss rate range 
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during the central He-burning phase, and our standard assumptions on the binary metalhcity. The 
models lose ~ 2Mq during the Hc-burning phase, and the mass loss cfRcicncy is so large that the 
stellar models show a very anomalous path in the H-R digram. Such a high mass-loss rate - if con- 
clusively established and, moreover, if typical of Cepheids also in the LMC - would be compatible 
with standard stellar evolution models only in case of episodic mass-loss events, that do not change 
the total mass of the star by more than a few hundredths of solar masses. 

3. Conclusions 

The recent study by Pietrzyhski et al. (2010) of a classical Cepheid in the well detached, double- 
lined, eclipsing binary OGLE-LMC-CEP0227, has provided the first very accurate determination 
of the dynamical mass of a classical Cepheid variable, and its radius. Our analysis has shown that 
stellar evolutionary models widely employed to study Galactic and extragalactic stellar populations 
are able to match - assuming the typical metallicity of LMC Cepheids - the derived mass and radius 
of both components with a single value for the age of the system, without any specific fine-tuning. 
Our result implies that there is no discrepancy between dynamical and evolutionary masses for the 
Cepheid star in this system, when current stellar models - calculated with up-to-date stellar physics 
and extensively tested against independent empirical benchmarks - are used in conjunction with 
the method employed in this analysis. 

It is also worth noticing that the agreement between evolutionary masses and high precision 
dynamical masses of Cepheids discussed in this paper, is based on just a single eclipsing binary 
system. The discovery and analysis of more systems similar to OGLE-LMC-CEP0227 - in con- 
junction with direct measurements of their chemical composition - will be enormously beneficial 
to the fields of stellar evolution and pulsation theory. 
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Table 1: Parameters of OGLE-LAIC-CEP0227 





Primary (A) 


Secondary (B) 


Mass (M/Mq) 


4.14±0.05 


4.14±0.07 


Radius (R/R©) 


32.4±1.5 


44.9±1.5 


Teff(K) 


5900±250 


5080ib270 
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solid - Mi= 4.14Mq, M,= 4.1 IMq I 
long dash - U.= 4.17Mg„ M,= 4.14Mq I 
short dash - M,= 4.21Uq, M,= 4.18Mq J 



I I I I I I I I I I I I I I I I 

145 150 155 160 

Time (Myr) 

Fig. 1. — The age - radius diagram: empirical values for the radii of the primary (labelled as A) 
and secondary (labelled as B) components are displayed as solid horizontal lines. The horizontal 
dash-dotted lines mark the boundaries of the la error regions. Three theoretical models with 
different values of the initial main sequence mass (Mj) and corresponding mass during the He- 
burning phase (Mf) are also displayed. The mass loss is efficient mainly along the red giant phase 
before central He-ignition, and reduces the stellar mass by a constant amount, equal to 0.03Mq. 
The three selected values of Mi and Mj are all within the la error bar of the empirical estimates 
(see Table 1). 
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Fig. 2. — The radius - effective temperature diagram: the stellar evolution models are the same as 
in Fig.l, while the filled circles and associated error bars correspond to the empirical estimates for 
the binary system components. The location of both the Fundamental Blue Edge and Fundamental 
Red Edge of the Cepheid instability strip are also plotted. 



